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Abstract 

a 

Using laser tweezers, we study colloidal interactions of solid microspheres in 
the nematic bulk caused by elastic distortions around the particles with strong 
tangential surface anchoring. The elastic interactions overcome the Brownian 
motion when the interparticle separation is less than 3 particle diameters. 
The particles attract when the angle 6 between and the uniform far-field 
director no is between 0° and ~70° and repel when 75° < 9 < 90°, which 
£T) \ is different from the model of quadrupolar interactions expected for large 

separations. The particles aggregate in chains directed at ~30° to no and, at 
■ higher concentrations, form complex kinetically trapped structures. 
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Interaction of particles in the bulk [1-5] and at the surfaces [6] of nematic liquid crystals 
(LCs) provide a rich variety of new physical phenomena of interest to both fundamental and 
applied science of liquid crystals and colloids. The studies recently expanded into the area of 
optical manipulation of LC colloids [7-9] and even to the development of biological sensors 
[10]. Anisotropy of molecular interactions at the particle surface leads to elastic distortions 
of the LC director n. A spherical particle immersed in LC with a uniform far-field director no 
creates distortions that are usually of dipolar or quadrupolar symmetry [1-6]. In the dipolar 
case, the interactions lead to chaining of particles along n [1,4,5,7-9,11] and a number of 
more complex structures [6,11-14]. Further progress is impossible without a quantitative 
characterization of the elasticity-mediated forces. However, the data are available only for 
the case of dipolar interactions of particles with normal surface anchoring of n [7,9,12]. 
Tangentially anchored spheres are of especial interest since they are expected to interact as 
quadrupoles [2,3]. Their interaction has never been characterized quantitatively, although 
it has been shown that water droplets attract each other along the direction 9 = 30° and 
eventually coalesce in the nematic LCs [11]. Here we report on the first direct measurement 
of colloidal interactions for spheres with tangential boundary conditions in the nematic LC. 
We employ fluorescence confocal polarizing microscopy (FCPM) [15] to characterize the 
director distortions around the particles and optical trapping with laser tweezers to measure 
the pair interaction force ~P P as a function of the separation between their centers and 
the angle 9 between and the far- field director n . 

A proper choice of materials is of prime importance for both FCPM and optical tweez- 
ers experiments. We use a low birefringence (An = 0.04) nematic ZLI2806 (EM Chem- 
icals) doped with 0.01wt.% of fluorescent dye n,n'-bis(2,5-di-tert-butylphenyl)-3,4,9,10- 
perylenedicarboximide (BTBP, Aldrich) for FCPM [15]. We use fluorescently-labeled 
Melamine Resin (MR, Aldrich) spheres of diameter D = 3um. The refractive index of MR 
n — 1.68 is higher than the average refractive index n = a/ (2tiq + n 2 ) /3 ~ 1.49 of ZLI2806 
(no is the ordinary and n e is the extraordinary refractive indices) to assure stable particle 
trapping [9]. The beads were coated with thin polyisoprene layers to produce tangentially 
degenerate alignment with a vanishing azimuthal anchoring coefficient W a J$ 10~ 10 J/m 2 
and polar anchoring coefficient W p = (0.6 ± 0.4) x 10 -4 J/m 2 [16]. The cells are formed by 
thin (0.15 mm) glass plates; the cell gap h = (30 — 100)um is set by mylar spacers. Thin 
polyimide layers (PI2555, HD MicroSystems) at the inner surfaces of the plates are rubbed 
to set a uniform in-plane director n . 

We use a fast version of FCPM based on a Nipkow-disk confocal system integrated with 
the Nikon microscope Eclipse E-600. We determine the director configurations (through 
polarization-dependent fluorescence of BTBP) as well as the positions of particles (through 
fluorescence of Rhodamine B labelling the MR spheres) [9] . We use an optical manipulator 
(Solar-TII, LM-2) and a laser (A = 1064nm, used power P = — 200mH^) for the dual- 
beam laser trapping [9]. The optical trap is formed by a lOOx objective (NA = 1.3) [17]. 
Quantitative measurements of forces in LCs require some additional care for three reasons: 
(1) the difference in refractive indices of the particle and the host depends on the local n, (2) 
the focused light beam can reorient n [8], and (3) light defocusing in a birefringent medium 
can widen the laser trap [9]. We mitigated all these problems by using a low-birefringent 
LC and colloidal particles larger than the waist of the laser beam 0.8/xm) [9]. 

Isolated spheres with tangential anchoring create distortions of the quadrupolar type 
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that quickly decay with distance r from the center of the particle and become optically 
undetectable at r ^ 2D, Fig. la. The director tangential to the sphere forms two point 
defects-boojums at the poles of particles [18]. When the beads are free to move around, they 
attract and form chains oriented at 9 xs ±(25° — 35°) with respect to n , Fig.lb,c. The chains 
elongate with time as other beads join. To verify whether the result is not an artefact of the 
possible tilt of the chain with respect to the cell plane, we used FCPM observations of the 
vertical cross- sect ions of thick samples with chains well separated from the bounding plates, 
Fig. Id. FCPM shows that the aggregation angle is indeed 9 ~ ±(25° — 35°), with most 
of the particles joined at 9 m ±30°. The same angle is observed not only in ZLI2806 with 
nearly equal splay K 1 and bend K 3 elastic constants, but also in pentylcyanobiphenyl (5CB) 
with K 3 /Ki w 1.56, in the mixtures E7 with K z jK x w 1.5 and ZLI3412 with K z jK x w 1.1; 
droplet attraction at 9 ps 30° has been also reported by Poulin and Weitz [11]. 

The chains in the bulk can glide on conical surfaces with the apex angle 29 ~ 60°, 
Fig.le,f. If the concentration of particles is sufficient, the chains interact and form more 
complex structures, Fig.le,f. Chain aggregation results in kinetic trapping (jamming) of 
particles, similar to the case of isotropic colloids [19,20] but different in the sense that the 
structures in LC reflect the anisotropy of interactions, Fig.le. The anisotropic features of 
the formed structures persist in a broad range of concentrations, Fig.lb,e. Unless the sample 
is heated to the isotropic phase, the jammed structures do not break apart for months. 

To characterize the pair interaction, we used a dual beam laser trap. Two particles are 
trapped at different locations and then released by switching off the laser. Their motion 
is strongly influenced by Brownian motion when r p is large (pair #1 in Fig. 2a), but the 
anisotropic colloidal interactions become noticeable at r p < (3 — A)D. The particles attract 
if r p makes the angle 0° < 9 < 70° with n and repel if 75° < 9 < 90° (pair "2" in 
Fig. 2a). The attracting particles eventually touch at 9 ps 30°; the particles in the 0-range 
of repulsion continue to undergo fluctuations until they enter the zone of attraction. We 
further probe the pair interactions by releasing only one bead from the trap while keeping 
the second trapped; the original separation is kept constant, r p = 7.5/xm, Fig.2b. The drift 
of the particles reflects the direction of interaction force. 

To get a better insight into the pattern of repulsive/attractive forces, we perform an 
experiment in which the position of one particle is fixed by a strong trap (P = lOOmW) 
and the second particle is slowly (~ fim/s) moved around it by a low-intensity trap (Pm < 
5mW), Fig. 3. The second trap follows a circle of radius r t centered at the first trap. The 
colloidal interactions make the actual separation r p of the particles different from r t : r p < r t 
when the particles attract and r p > r t when they repel. For relatively large distances, 
r t > 2D, the strongest attraction is observed at 9 — 35° —45°, Fig.2a,3, which is higher than 
the aggregation angle 9 ~ 30°. 

To determine F p we first calibrate stiffness a of the optical trap for different laser powers 
P [21]. We measure the displacement Al d of the trapped bead from its equilibrium position 
in response to the viscous drag force F d = 3nDr]^ e ffV as the particle is moved at a constant 
speed \? || n in the LC with effective viscosity r)\\ e ff = 5.9- 10~ 5 Pa -s [22]; the trap stiffness 
is a(P) = F d /Al d (P). The amplitude of ~P P is then determined by the displacement Al p 
caused by the pair interaction, F p = aAl p [21]. We also use the technique of particle escape 
to probe interactions at r p > 1.5D [7,9]. A single particle is moved by the tweezers with 
an increasing velocity V || n until it escapes from the trap at which moment the "trap 
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escape" force is determined by the Stoke's law, F te = 37rDr]\\ e ffV. F te is found to be a linear 
function of P, indicating that nonlinear light-induced processes (such as realignment of n 
by the laser beam [8]) are insignificant [17]. F p is measured as F te when the escape from 
the trap is caused by the pair interaction [7,9]. The results obtained by the two methods 
above are in a good agreement. The advantage of both approaches is that the effects of 
anisotropic viscosity of LC are avoided (the velocities of the beads during measurement of 
F p are negligibly small) and that both a and F te are practically independent of n around 
the beads (An is small). 

To collect the data in Fig.4, the particles are placed at positions with different #'s but 
constant r p , and the force is determined by measuring Al p . When the laser power P is 
reduced, the direction of bead displacement (determined from the optical microscope images) 
indicates the direction of F p . Figure 4b shows how the angle 7 between ~r* p and r p changes 
when one particle circumnavigates the other. We measured the distance dependence of 
Fp( r p) f° r 9 — 30° = const, Fig. 5. At r p <J 1.5D, the dependence F p (r p ) is close to a power 
law oc r~ 6 , i.e., it is much sharper than in the case of spheres with dipolar distortions of n 
[5,12]. The force amplitude is of the order of K for r p < 2D (K « !5pN for ZLI2806). 
In addition to the experiments with D = 3/im, we measured the forces for larger particles, 
D = (4 — 7.5) urn. The qualitative picture remains the same, but Fp increases with D; for 
example, ~F* P increases by ~ 40% when D increases from 3um to 4um. 

Let us compare the experimental ~P p (6,r p ) to available theoretical models for quadrupo- 
lar interactions with the pair potential derived in the approximation K\ = K 3 = K [2]: 

U q (r P , 9) = ^ _ _E_ j (9 + 20 cos 29 + 35 cos 40) , (1) 

where R = D/2. In Fig.4 we plot the theoretical predictions for J? p (9,r p = 1.5D) and the 
^-dependence of the angle 7 between and ~P P = —VU q as calculated using Eq.(l) with 
K = 15pN and W p = 5 • 10~ 5 J/m 2 ; the calculated force is of the same order as in the 
experiment. The model predicts maximum attraction at 9 ~ 49° which is independent of 
r p , while the experiment shows maximum attraction at 9 ~ 30° for D < r p < 1.7D and 
at 9 pa 35° — 45° for 2D < r p < AD. Furthermore, the experimental repulsive sector is 
limited to 75° < 9 < 90°, while the model predicts repulsion also for 9^0. The model 
dependence F p oc r p 6 holds only for r p > 1.5D; at smaller distances it is much weaker, Fig.5. 
The discrepancies are not surprising, as the model [2,3] assumes r p » D and weak surface 
anchoring, K/W p » D. In the experiment, the elasticity-mediated forces overcome the 
background Brownian motion only when r p < (3 — 4)D; besides, K/W p ~ 0.2/im, smaller 
than D. Note that as r p /D increases, the agreement between the experiment and the 
model [2,3] improves, Figs. 2-5. To describe the interactions at close separation, the theory 
should take into account finite surface anchoring and deviations of the director from the 
quadrupolar to lower symmetries, similar to the case revealed in simulations of antiparallel 
elastic dipoles [4] . One possible reason for the lower symmetry of director distortions might 
be non-sphericity of the particle shape and/or non- homogeneity of the surface anchoring, 
even when the particle is perfectly spherical, see, for example Ref. [24]. 

To conclude, we characterized the pair interaction of colloidal spheres in the nematic 
host that induce quadrupolar director distortions when isolated. The interactions deviate 
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from the quadrupolar model because the experimental forces are significant only for small 
inter-particle distances. When the concentration of particles increases, they aggregate, first 
in chains tilted with respect to the far-field director and then into more complex aggregates. 
The chains break apart once the LC is heated to the isotropic phase (provided that the 
beads are charge-stabilized); the system, therefore, is of interest to study kinetic arrest and 
jamming [19] in anisotropic media with both attractive and repulsive interactions [20]. 

We thank B. Senyuk, S. Shiyanovskii, and H. Stark for discussions. We acknowl- 
edge support of the NSF Grant DMR-0315523 (Kent) and the AFOSR DURINT grant 
F496200110358 (SUNY). U.S. acknowledges the support of the Fellowship of the Institute 
of Complex and Adaptive Matter. 

FIG.l. Colloidal aggregation in the nematic bulk: (a) FCPM texture of the director 
distortions around a pair of particles; (b) FCPM image of director distortions in the LC co- 
localized with the fluorescence signal from particles (green) forming chains at 9 w 30° to n ; 
(c) scheme of the director distortions around the spheres in chains; (d) FCPM vertical cross- 
sections along the d-d line in (b) illustrating that the particle chains (green, fluorescence 
from Rhodamine B) are located in the bulk of LC (gray-scale, BTBP); (e) aggregation 
of spheres leading to kinetically-trapped network that reflects anisotropy of the host even 
at high particle concentration m 5 wt.%; (f) scheme of chain orient at at ions. The FCPM 
polarization is perpendicular to n in (a,b) and along the chain in (d). 

FIG. 2. Trajectories of particle pairs in the LC: (a) both particles are released from the 
optical traps at different locations (labelled by pairs of numbers 1 through 8; the trajectories 
of the two particles released at the same time are marked by the same color); (b) one particle 
is fixed and the other is released from a semicircle of radius r p pa 7.5/j.m at different 6. 

Fig. 3. Angular dependence of the center-to-center distance r p between two particles 
measured when one trap moves around the other trap along a circle of a radius r t . 

FIG. 4. Angular dependencies of for two MR particles in the LC: (a) F p {0) measured 
for r p = 1.5D (triangles), r p = 1.75D (squares), and r p = 2D (circles); (b) the angle 7 
between F^p(#) and as a function of 6; (c) vector representation of as determined 
experimentally for r p = 1.5D and r p = 2.5D and calculated for r p ^> D using (1). 

FIG.5. F p vs. r p measured for 9 « 30°. The solid lines are plots of F p = U8(D/r p ) 6 pN. 
The inset is the log-log plot of the same data. 
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